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In the absence of phonon contribution, a weakly coupled single orbital non-interacting quantum 
dot thermoelectric set up is known to operate reversibly as a Carnot engine. This reversible op- 
eration, however, occurs only in the ideal case of vanishing coupling to the contacts, wherein the 
transmission function is delta-shaped, and under open-circuit conditions, where no electrical power 
is extracted. In this paper, we delve into the thermoelectric performance of quantum dot systems 
by analyzing the power output and efficiency directly evaluated from the non-equilibrium electric 
and energy currents across them. In the case of interacting quantum dots, the non-equilibrium 
currents in the limit of weak coupling to the contacts are evaluated using the Pauli master equation 
approach. The following fundamental aspects of the thermoelectric operation of a quantum dot set 
up are discussed in detail: a) With a finite coupling to the contacts, a thermoelectric set up always 
operates irreversibly under open-circuit conditions, with a zero efficiency, b) Operation at a peak 
efficiency close to the Carnot value is possible under a finite power operation. In the non-interacting 
single orbital case, the peak efficiency approaches the Carnot value as the coupfing to the contacts 
becomes smaller. In the interacting case, this trend depends non-trivially on the interaction pa- 
rameter U. c) The evaluated trends of the maximum efficiency derived from the non-equilibrium 
currents deviate considerably from the conventional figure of merit zT based results. Finally, we 
also analyze the interacting quantum dot set up for thermoelectric operation at maximum power 
output. 

PACS numbers: 73.63.Kv,85.35.Gv,85.80.Fi,84.60.Rb 



I. INTRODUCTION 

Thermoelectrics are currently an object of immense in- 
terest and intense research activity owing to the possible 
enhancement of the energy conversion cfRcieiicy via nano- 
structuring [H-Q and novel materials design A higher 
energy conversion efficiency of a thermoelectric system is 
typically quantified by an increase in zT, the dimension- 
less figure of merit. The figure of merit zT is defined 
as: 



zT = 



(1) 



where S,a, and Kei^(ph) are the linear response trans- 
port coefficients, namely the thermopower (Seebeck co- 
efficient), the electrical conductivity and the electron 
(phonon) thermal conductivity, with T being the average 
operating temperature. The proposed increase in zT is 
envisioned via novel approaches towards engineering the 
electronic [H-Q or phononic transport 043 properties. 
Among various low-dimensional nanoscale systems, zero- 
dimensional systems such as quantum dots have been of 
special interest because they may exhibit an infinitely 
high value of zT in the absence of phonon thermal con- 
ductivity 

The energy conversion efficiency 77 of a thermoelectric 
system is usually defined as 77 = with P being the 

extracted power and J'q being the input heat current. 
Consider, for example, a set-up with a central system 
sandwiched between two reservoirs held at a fixed tem- 
perature and electrochemical potential. Under the as- 
sumption of small electrochemical potential, A/i, and 



small temperature, AT, differences between the reser- 
voirs, the electric current J, and the heat current Jq 
may be written as H, : 



J^LiiA^j + LisAT, 

Jq^ L21AH + L22AT, 



(2) 



where Lij represent the Onsager coefficients. The On- 
sager coefficients Lij, are in turn related to the linear 
response parameters, namely cr, S, and Kg/, that appear 
in the aforementioned definition of zT. The efficiency 77, 
when maximized with respect to J, yields its maximum 
Vmax to be an increasing function of zT [Tl| given by: 



Vc 



J H 



(3) 



where the thermoelectric material operates between two 
contacts maintained at temperatures Th and Tc, with 
77c = 1 — ^ being the Carnot efficiency. It is therefore 
convenient to employ zT as a performance metric to fa- 
cilitate the design of maximally- efficient thermoelectrics. 

However, as noted above, the use of zT as the per- 
formance metric in lieu of the actual efficiency relies on 
the assumption of linear response. While a high figure of 
merit zT is often a necessary component for a good ther- 
moelectric, it docs not sufficiently underpin the working 
conditions that are involved. For example, an analysis 
of the figure of merit zT of the single orbital quantum 
dot system Q , under the condition of vanishing coupling 
to the contacts, simply points to its infinite value and 
the resulting efficiency maximum as the Carnot value. It 
was pointed out only recently p^ - [l7j that this efficiency 
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maximum only occurs under open circuit conditions, im- 
plying an operating condition with a vanishing current 
and hence a vanishing power output. The Carnot effi- 
ciency is reached only due to the possibility of this re- 
versible operation [l2,[l4| under open circuit conditions. 

The open circuit condition, although associated with a 
vanishing current, is an operating point which has both 
an electrochemical potential gradient and a temperature 
gradient. The voltage Vs associated with this electro- 
chemical potential difference ±qVs, with q being the elec- 
tric charge, is known as the Seebeck voltage. This volt- 
age cancels the current set up by the applied temperature 
gradient. It is hence pertinent to analyze thermoelectric 
operation by using a nanocaloritronic set up, wherein the 
central system is subject to a bias drop, not necessarily 
equal to the Seebeck voltage, and a temperature gradi- 
ent. 

The central system considered in our nanocaloritronic 
analysis is a single orbital interacting quantum dot. 
Steady state non-equilibrium currents through the cen- 
tral system, rather than linear response parameters, are 
used to evaluate the power and hence the efficiency at 
each operating point. Each operating point is defined 
by the applied bias and the applied temperature gra- 
dient. The primary goal of our transport calculations 
is to identify the operating conditions that point to a 
specific operating efficiency in relation to the operating 
power. Some recent works [ij, ElHIil, for example, 
have specifically analyzed the operation of a single or- 
bital non-interacting quantum dot thermoelectric set up 
at maximum power. A recent investigation that includes 
Coulomb interactions [l^ has noted the importance of 
non-linear effects, and has focused on the role of a phonon 
bath on the thermoelectric operation. Another recent 
work [2^ has focused on the effect of Coulomb interaction 
on the figure of merit zT. The main focus of this paper 
therefore is a comprehensive performance analysis of a 
quantum dot thermoelectric set up. The following funda- 
mental aspects of the thermoelectric operation of a quan- 
tum dot set up arc then discussed in detail: a) With a 
finite coupling to the contacts, a thermoelectric set up al- 
ways operates irreversibly under open-circuit conditions, 
with a zero efficiency, b) Operation at a peak efficiency 
close to the Carnot value is possible under a finite power 
operation. In the non-interacting single orbital case, the 
peak efficiency approaches the Carnot value as the cou- 
pling to the contacts becomes smaller. In the interacting 
case, this trend depends non-trivially on the interaction 
parameter U. c) The evaluated trends of the maximum 
efficiency derived from the non-equilibrium currents de- 
viate considerably from the conventional figure of merit 
zT based result. Wc point out in detail the discrepan- 
cies between our non-equilibrium analysis, and the linear 
response analysis that is usually based on the figure of 
merit zT . Civen the current experimental possibility of 
thermoelectrics across zero-dimensional systems [2l|, , 
and the recent theoretical activity exploring non- linear 
thermoelectric effects [H, [H, [H, [l^, [3, [3 across them. 



our paper elucidates the importance of Coulomb interac- 
tion on their thermoelectric performance. 

This paper is organized as follows. Section II describes 
the necessary formulation: first the definition of the elec- 
tric and energy currents through the quantum system, 
then the formalism used to evaluate these currents and 
hence the power output and efficiency across it. The 
quantum transport system under consideration is a single 
level Anderson-impurity type quantum dot that is weakly 
coupled to the contacts in the sequential tunneling limit. 
The formulation for currents follows from the density ma- 
trix master equation a ppr oach under this sequential tun- 
neling approximation [23|-|26|. Section III begins by de- 
scribing the thermoelectric operation of a quantum dot 
set up in the absence of interactions ([/ = 0). Follow- 
ing that, the fundamental results due to the introduction 
of Coulomb interactions (finite U) are discussed in de- 
tail. The section concludes with an analysis of the max- 
imum power operation. It is shown that with Coulomb 
interactions the maximum power operation is relatively 
unaffected in comparison with the non-interacting case 
discussed in other works [l^ [l^, [l^ . Section IV summa- 
rizes the results of this work. 



II. THEORETICAL FORMULATION 

A prototype nanocaloritronic configuration of a quan- 
tum thermoelectric set up is shown in Fig. [Ija). In this 
set up, the central quantum system is coupled to two 
reservoirs/contacts a, which are labeled a = H, (hot) 
and a = C (cold), each of which is characterized by 
a temperature Tfj{c) and an electrochemical potential 
f-H(c) ■ This set up closely relates to that of a heat engine 
commonly studied in classical thermodynamics. Under 
this set up, the central system plays the role of the heat 
engine with the reservoirs acting as heat sources or sinks. 
In the case of a thermoelectric set up, however, we also 
need to invoke the additional concept of particle exchange 
[1^ , due to the fact that the reservoirs are characterized 
by both temperature and electrochemical potential, thus 
facilitating both energy and particle transport. An ap- 
plied voltage bias qVapp = fJ-c ~ fJ-H, an applied temper- 
ature gradient AT = Th — 7c . or both, trigger particle 
flux and hence a flow of both an electric and an energy 
current results. In order to describe quantum thermo- 
electric transport across the system of interest, electric, 
and energy currents must be clearly deflned. 



A. Electric and energy currents 

We begin with the fundamental thermodynamic equa- 
tion relating the internal energy E of a non-magnetic sys- 
tem with extensive variables such as entropy S, volume 
V, and particle number N given by |27| : 



E{S, V, N) ^TS ~pV + fiN, 



(4) 



where the intensive variables are the temperature T, pres- 
sure p, and the electrochemical potential /x. The above 
definition relates to the thermodynamic state of the sys- 
tem. The case of thermoelectric transport involves elec- 
tron transfer processes during which the thermodynamic 
state of the system changes. When such processes are 
involved, one measures the change in the internal energy 
with respect to the change in the extensive parameters. 
Specific to our case, thermoelectric transport occurs at 
constant volume {AV = 0). Electron transfer processes 
occur between either reservoir (a = H,C), each main- 
tained at a fixed temperature and a fixed electrochemical 
potential, and the system as shown in FiglTJa). One may 
then write an equation for the infinitesimal change in the 
internal energy due to an infinitesimal charge transfer 
between either reservoir (a = H, C) and the system as: 



dEa = TadSa + IJI-adNa 



(5) 



a) 



b) 



fir kj 



2s + U 
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P^-JV 



c) 



fir kj 



H C 



Using the above definition, we can take total time deriva- 
tives to define a current associated with the correspond- 
ing flux of the extensive variables given by: 

~ TaJg -|- flaJ]\f (6) 

The quantity TaJg is usually termed as the heat cur- 
rent involved in the isothermal electron transfer between 
either reservoir and the system. This quantity, in gen- 
eral, is the contribution to the energy current that keeps 
track of entropy flow as given by Js- Although the term 
heat current is widely used in literature, it may not by 
itself be accurate because heat is not a state function 
and the definition of differentials of such quantities may 
not be obvious. However, in the case of thermoelectric 
transport that is considered in this work, the quantity 
Jq may be termed heat current, fbllowing the arguments 
put forward in early works hOl. The central assump- 
tion is that the reservoirs are maintained in equilibrium 
and hence the flow of charge and heat to and from the 
reservoirs happens reversibly. All irreversible processes 
are expected to occur in the interfacial region between 
the reservoirs and the quantum dot. In such a reversible 
process, the expression AQa = TaASa holds true, for 
the reservoir. Hence, the equation may now hence be 
recast in terms of the quantity Jq = TaJg as: 

where, the quantity Jq is the heat current in the reservoir 
a. In the quantum mechanical case, in order to deflnc 
currents, we define the time-dependent average current 
due to an operator O that is associated with one of the 
extensive variables as follows: 



Jo{t) = 




FIG. 1: (Color online) Nanocaloritronics of a quantum ther- 
moelectric transport set up. a) A typical thermoelectric set 
up comprises of the central quantum system described by the 
Hamiltonian Hs sandwiched between two reservoirs labeled 
hot(cold) a = H{C). When this central system is subject 
to an electrochemical potential gradient and a temperature 
gradient, the resulting current J drives an electrical power 
P = JVapp via the electrical leads. Equal contact couplings 
Th = Fc are assumed throughout, b) A single orbital quan- 
tum dot is parameterized by its single particle energy level e 
and the Coulomb interaction parameter U . Transport is rep- 
resented as transitions between states of the many-particle 
spectrum with electron numbers differing by ±1. Transport 
channels then comprise the energy difference e,e -\- U be- 
tween those states with electron numbers differing by ±1. c) 
Schematic depicting the thermoelectric effect under open cir- 
cuit conditions: The built-in or Seebeck voltage Vs enforces 
zero current in the circuit. Useful power can be extracted 
when the applied voltage Vapp G [0, Vs] , where the condition 
Vapp = Vs enforces open circuit operation. The thermoelec- 
tric efficiency defined in the operating region Vapp £ [0, Vs] is 
strongly affected by the energy difference e — fia, the applied 
temperature gradient AT = Th — Tc, and the magnitude U 
of the Coulomb interaction. 



where [H ^ O] represents the commutator of the overall 
Hamiltonian H with the operator O. For a quantum 
mechanical set up, based on the schematic of Fig. [TJa), 
the description of currents thus begins by describing 
the overall Hamiltonian H which is usually written as 
H = Hs + Hr + Hti where Hs, Hn and Ht represent the 
system, reservoir and reservoir-system coupling Hamilto- 
nians respectively. In this paper, the system comprises of 
the single orbital Anderson impurity-type quantum dot 
subject to Coulomb interaction described by the follow- 
ing one-site Hubbard Hamiltonian: 
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where represents the orbital energy, ha = d^^da is 
the occupation number operator of an electron with spin 
a =1, or a and U is the Coulomb interaction be- 
tween electrons of opposite spins occupying the same or- 
bital. The exact-diagonalization of the system Hamil- 
tonian then results in four Fock-space energy levels la- 
beled by their total energies 0, e-j-, and e-|- -I- ej. -I- \J . 
In this paper, we consider only a spin-degenerate level 
such that e = e-f = e|.. Electronic transport gener- 
ally involves the addition and removal of electrons. In 
the limit of weak contact coupling (^iF ^ fc^T), trans- 
port may be viewed as of transitions between the Fock- 
space levels that differ by an electron number of ±1 
as shown in Fig. [ijc). The reservoir/contact Hamil- 
tonian is given by = Y.a=H,cY.ka '^o.kahaka = 

Tliaka ^aka, whcrc a labcls the hot/cold reservoir {H or 
C in our case) and the summation is taken over the sin- 
gle particle states labeled {ka}. The tunneling Hamil- 
tonian represents the system-contact coupling usually 

written as Ht = Eafc<T {tckaC^f^Ja + t*^kJt^c. 

Tliaka^Taka, whcrc (ct,c) and {d\d) denote the cre- 
ation/annihilation operators of the reservoir and system 
states respectively. 

Pertinent to our problem, one can use Eq.® to eval- 
uate, for example the steady-state electric, and energy 
currents through the system. The steady state current 
is then derived in the limit when t — > oo. Also, in our 
case the operator does not explicitly depend on time im- 
plying that ^ = 0. We can then write the particle 
current due to either contact a = H/C by summing 
contributions over its one-electron states labeled {ka} 



Yka 



as J' 



N 



ka dt 



■). The expression 



for the electric current due to either contact given by 
J" = ^iJn then becomes: 



time derivative of the number operator haka involves the 
evaluation of its commutator with the tunneling Hamil- 
tonian hTaka, as described in Ea.([TUl). Following Eq.®, 
the expectation value of an operator is evaluated by trac- 
ing over the composite system-reservoir density matrix 

i.e., {^O^ = TT{p{t)d{t)}. The time evolution of p{t) is 

given by the Liouville equation. The reduced density ma- 
trix Pred(t) of the system may be obtained by performing 
a trace exclusively over the reservoir space. An expan- 
sion of the Liouville equation to the second order in the 
tunneling Hamiltonian in the limit of weak contact cou- 
pling {hT ksT), leads to the density matrix master 
equation for the reduced density matrix of the system 
[23l - l25l |. In second order, coherences vanish for the con- 
sidered single orbital model [1^, and one is left with a 
scalar rate equation [28l - l30j in terms of the occupation 
probabilities Pf^ = {N,i \ Pred{t) \ N,i) of each TV elec- 
tron Fock state \N,i) with total energy Ej^ . The index 
i here labcls the states within the iV electron subspace. 
This Pauli-mastcr equation then involves transition rates 
R{N,i}^{N±i.j} between states \N,i), and jA^ ± l,j) dif- 
fering by a single electron, leading to a set of independent 
equations defined by the size of the Fock space: 



dP, 



N 



dt 



N±lJ)^{NA)Pj 



N±l 



R 



(13) 

along with the normalization equation J2i n -^1^ ~ ^■ 
Notice that, in the stationary limit considered here, 
where < — >■ oo, the Markov approximation implicit in 
Eq. (fT3|) becomes exact [1^ HB] • At energies close to the 
Fermi level, metallic contacts can be described using a 
constant density of states, parameterized using the bare- 
electron tunneling rates 7^ = J2ka ^\taka,s\'^S{E - Cka), 
with {a = H/C). Wc define the rate constants as: 



J"' = ^^('^'^--Z [hTak'a' ,hakcr] 



h(7 h' \ 



(10) 



where q is the electronic charge. Likewise, the energy 
current due to the contact a is written as: 



dt 



E E ^Ti^Tak'a' ,hakcr] 



ka k'a' 



E ^otka- 



dfl 



aka 



ka 



dt 



(11) 



The quantity Jq due to the contact a then becomes: 



Jq — TaJg — ( ^ ^ (EgfctT ~ Ma) 



d^aka 

dt 



(12) 



The above expression is the commonly employed relation- 
ship that connects the so called heat currents with en- 
tropy and particle currents [l5| . Calculation of the total 



^^\{N,i\d„\N + l,3)?. 



aij 



(14) 



The transition rates for the removal (|A^, i) |A^— l,j)), 
and addition {\N,i) — > \N + 1, j)) transitions arc then 
given by 



R 



(7V,i)^(JV-lj) 



E 



Nr 
ij 



a=H,C 



l-f 



Nr _ 



knT„ 



R 



^ aij J 



fJ-a 



a=H,C 



ksTa 



(15) 



The contact electrochemical potentials and temperatures 
are respectively labeled as pa and Tq, and / is the cor- 
responding Fermi-Dirac distribution function with single 
particle removal and addition transport channels given 
by 



Nr 
Na 



E. 



N 



E 



N-l 



Ef-^ 



E. 



N 



(16) 
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Finally, the steady-state solution to Eq. p3|) . set by 

= 0, is used to obtain the terminal current asso- 
ciated with contact a: 

Ntat 

N=l ij 



-Hr, 



(17) 



where Ntot is the total number of electrons in the system. 
In our case, for example, Ntot — 2. Likewise, the quan- 
tity Jq, associated with either contact can be similarly 
defined using Eq. p2)) as: 



Ntot 



N-1 



N=l ij 



Here, the sum over reservoir indices (ka) in Eq. p^ . has 
been replaced by indices («, j) corresponding to the sys- 
tem states because of elastic electron transfer between 
the reservoir and the system, described by the energetics 
fifecr = e-jf for the additive transition, and eka- = , 
for the removal transition. Notice from Eq. fTS)) . that the 
total rates R{N,i)^(N±i,j)i and R{N±i,j)^(N,i) appear- 
ing in Eq. (jl3p , are the sum of individual rates associated 
with either contact in Eq. (fT7| and (fTS]) . 



B. Power and efficiency 

In a classical heat engine, the efficiency of a thermo- 
dynamic cycle is defined as 77 = which is simply the 
ratio between the work extracted and the heat supplied. 
However, while working with the nanocaloritronic con- 
figuration shown in Fig. (Ha), it is important to evaluate 
the efficiency under a finite power operation because con- 
version of entropy currents to electric currents is desired. 
In our case, in which the operation at a finite power is 
desired, the efficiency is given in terms of the rates of 
flow of various quantities: 



Tin ' 



(19) 



where, the instantaneous power or just the power is de- 
fined as P = ( J^" - J§"*). Following Eq.H^, and assum- 
ing no intra-system or endo-dynamic energy changes due 
to inelastic processes, the net electrical power between 
the hot and cold reservoirs can be written as: 



P 



f^H)J = -VappJ, (20) 



where J = = — J*-^ refers to the electric current 
whose magnitude is conserved in steady state. It must be 
noted that the above expression has both the Joule (irre- 
versible) and the thermoelectric (reversible) components 



[3ll |33 | . For example, specific to the linear response case, 
one obtains P = Ln ((A//)^ -I- Li2ATA/i) /g, by em- 
ploying Eq. ([2|) . This combines linear and quadratic terms 
in the applied voltage bias qVapp = {^ic — Mi?) = A/x, 
the linear term being the thermoelectric part, and the 
quadratic term being the Joule part. 

The power generated, and hence the efficiency, is gen- 
erally evaluated at an operating point. Each operating 
point is specified by the applied bias Vapp and the tem- 
perature gradient AT = Th — Tq- For the upcoming 
analysis, we work with the convention that the temper- 
ature gradient is applied at the contact labeled H , and 
the voltage bias Vapp is applied at the contact labeled C. 
In all our calculations, we assume that half of the applied 
voltage drops across the quantum dot as a result of equal 
capacitive coupling to the two contacts. 



III. THERMOELECTRIC OPERATION OF A 
QUANTUM DOT 

In the realm of molecular electronics or quantum dot 
transport, it is common to start with a microscopic un- 
derstanding of transport processes across a single spin 
degenerate orbital subject to Coulomb interactions. Of- 
ten this leads to a qualitative physical picture of various 
experimental observations and the additional complexity 
of multiple levels may append mainly to the quantitative 
aspect. Based on the formulation discussed in the previ- 
ous section, we first elucidate the thermoelectric opera- 
tion of the quantum dot set up without Coulomb inter- 
actions. Following that, we discuss the important results 
of this work that arise due to the inclusion of Coulomb 
interactions. 



A. Power and efficiency of a non-interacting 
quantum dot thermoelectric set up 

First we discuss the results that follow from the se- 
quential tunneling model. This model implies a delta line 
shape for the quantum dot density of states and trans- 
mission function. In this limit, the analytical result for 
the currents are given by: 

J = ^^^^(,fe(e)-/c(e)). 



IH + IC 
IH + IC 



(e-;U„)(/H(e)-/c(e)), (21) 



with jH^c being the contact coupling energies associated 
with contacts H, C. Here, /a(e) = / ^f^f^) refers to the 
Fermi-Dirac distribution of either contact. The factor of 
2 appears due to spin-degeneracy in the non-interacting 
case. 

Based on the schematic in Fig. [IJc), the basic thermo- 
electric operation can be described as follows. An electric 
current is set up by the applied temperature gradient. 
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Under open circuit conditions, the Seebeck voltage Vs 
is set up in order to oppose this current. This built-in 
voltage can be used to drive power across an electrical 
system, say a resistor. Alternatively, while working in 
a circuit configuration an externally applied voltage bias 
Vapp may be used as a variable electric current source. 
The condition that enforces zero electric current is then 
equivalent to an operating condition with an applied bias 
Vapp = Vs. It is now easy to see from Ea.([2T|) that a zero 
electric current may be enforced by ///(e) = /c(e) or 



fJ-H 



ksTc 



(22) 



In general, the quantity Jg under the above condition 
need not also be identically zero. In the present case, 
however, Jq is also zero, and is easily noted from Eq. ([21]) . 
This point has an important implication with respect to 
the operating efficiency. 

It can be shown by using the definition of efficiency 
Ea. ([T^ and Eg. ([?!]) ■ that the efficiency is given by: 



ipc - Mg) 
(e-/-tff) 



(23) 



Notice that the expression for the efficiency is indepen- 
dent of the current J. Under open circuit conditions, 
it can then be deduced that the maximum efficiency 
Vmax = VC- Thus, a non-interacting quantum dot, under 
the limit of vanishing coupling to the contacts, operates 
reversibly and achieves the Carnot efficiency under open 
circuit conditions. A brief discussion of the thermody- 
namic aspects of this reversible operation is carried out 
in Appendix I. Finally, making another connection with 
[1], the zero value of the quantity Jq results in a zero 

electron thermal conductivity, i.e., n^i = = 0. This 
results in an infinitely high value of zT in the absence 
of phonon contribution, and following Eq.(l3]), leads to 
the maximum efficiency equaling that of the Carnot effi- 
ciency. 

Effect of line width: The rather surprising result of 
achieving a finite efficiency under zero power operation 
is indeed an artifact of the sequential tunneling approx- 
imation which implies the idealized delta form for the 
transmission function. Going beyond the sequential tun- 
neling approximation, the delta function broadens, and 
hence impacts the conclusions drawn above. In the spe- 
cific case of the non- interacting limit, it is possible to ex- 
actly evaluate the currents using for example, the trans- 
mission formalism Il3i as: 



J 



Jn — 



IH + 7C J- 



IH + IC 



/oo 
dED{E) {fH{E)~fc{E)), 
-OO 

/OO 
dED{E){E - ifniE) - fc{E)) , 
-OO 

(24) 
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FIG. 2: (Color online) Power and efficiency in the non- 
interacting {U = 0) limit for Carnot efficiency rjc = 0.23 
(black solid) and 77c ~ 0.33 (green circles). The tempera- 
ture at the cold contact is set to Tc ~ lOOK, and the equi- 
librium energy level placement is set to e — fiH ~ 2kBTH 
at Vapp ~ 0. The couplings to the reservoirs are taken as 
ft7i/ = h-yc = O.OlmeV ^ lO'^ksT. a) Plot of extracted 
power as a function of tlie applied bias Vapp- The span of the 
operating region Vapp G [0, Vs] broadens with increase in the 
applied temperature gradient AT — Th — Tc- Results from 
the sequential tunneling approximation (dotted) and the ex- 
act calculation (bold) are identical, b) Corresponding plots of 
efficiency in the operating region. Under the sequential tun- 
neling approximation (dotted), the efficiency maximizes at the 
Carnot efficiency rjc when the applied bias equals the built- 
in voltage [Vapp = Vs)- This corresponds to the reversible 
thermoelectric configuration [^. [T^. (see text). In the ex- 
act calculation, however (bold), the efficiency drops to zero 
under open circuit conditions. The efficiency at maximum 
power lies in an intermediate operating point corresponding 
to the maximum power Pmax shown in (a). 



where the broadened density of states D{E) is given by: 
1 



D{E) 



IH + IC 



e)' + ((7H+7c)/2)' 



(25) 



The efficiency, as seen in Fig. [^b), given by the ratio of 
P and Jq , drops to zero under open circuit conditions. 
This is because unlike in the previous case, it can be 
noted from Ea.([M)). that Jq need not also vanish when 
J = 0. This also implies that the Carnot efficiency can 
never be reached in the real situation. In our simulations, 
we have used a h'-^H = ^"1c — O.OlmeU = lO^'^fc^r, 
such that the condition for weak coupling to the con- 
tacts hV <^ ksT, is satisfied. One can hence note from 
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Fig. ^a) that there is almost no difference in the vari- 
ation of power between sequential tunneling approxima- 
tion and the exact calculation. The sequential tunneling 
limit thus provides a very good approximation for the 
evaluation of currents in the limit of weak coupling to 
the contacts, but fails to describe the correct trend for 
the efficiency in this limiting case. 

The efficiency under open circuit conditions is identi- 
cally zero, because the quantity Jq is finite when the 
current J is zero. Using Eq.([21]), and Jq — T^Jg, we 
note that: 

(26) 

where, j{E) = ^1M^D{E) {Jh{E) ~ fc{E)). Physi- 
cally, the above result implies that although the flow of 
electrons from the hot to the cold contact under open 
circuit conditions is balanced by the reverse flow, the net 
flow of entropy is not. Entropy can then be produced 
while maintaining a zero net particle flux. This entropy 
production thus results in a finite entropy current under 
open circuit conditions. Therefore, it implies that un- 
like the special case of vanishing coupling to the contacts 
that is discussed in Appendix I, spontaneous electron ex- 
change between the reservoirs is inherently irreversible. 
The entropy generated is dissipated as heat deep in the 
reservoirs. 

Operating region: Thus far, we have considered only 
one operating condition, namely that of the open cir- 
cuit operation in which Vapp = Vs- In order to fully 
characterize the thermoelectric set up, an understanding 
of its operation at an arbitrary applied bias Vapp must 
be considered. According to our convention, and follow- 
ing the definition of efficiency in Ea. ([T5|) . useful work 
may be extracted only in the region of positive power 
[P > 0). Therefore, the domain in which P > 0, that 
is represented by the applied bias < Vapp < Vs, as 
shown in the schematic in Fig. [IJa) , defines the operat- 
ing region. The extracted power in this operating region 
Vapp S [0, Vs] is shown in Fig.[2Ia), for two different val- 
ues of the applied temperature gradient, and hence of 
the Carnot efficiency. Notice that the extracted power is 
identically zero (P = 0) under two operating conditions: 
short circuit condition- when Vapp = 0, and open circuit 
condition- when Vapp = Vs ■ The first one corresponds to 
a zero bias and the second one corresponds to the con- 
dition with a zero electric current. The operating region 
also becomes larger as the applied temperature gradient 
AT, and consequently the associated Carnot efficiency 
rjc = is increased. This is because an increase in AT 
increases the amount of current flowing through the level, 
as a result of which a higher applied voltage Vapp = Vs 
is needed to counter it. The variation of power in the 
operating region is quasi-quadratic and has a maximum 
(marked Pmax) in the operating region. 

The efficiency in the sequential tunneling case increases 
monotonically and quasi-linearly (TtI . [isj from to rjc in 



the operating region. In the exact calculation, however, 
the efficiency reaches a maximum that is close to the 
Carnot value and then drops to zero at the open circuit 
operation point Vapp = Vs- The abruptness of this be- 
havior depends on how large the coupling to the contacts 
is. This deviation of the maximum efficiency obtained 
via the exact calculation from the ideal Carnot value ob- 
tained via the sequential tunneling approximation will 
become more pronounced as the contact coupling is in- 
creased. 



B. Power and efficiency of an interacting quantum 
dot thermoelectric set up 

With the same initial configuration as in the previous 
case, we now study the effect of varying U. Referring to 
the state transition diagram in Fig. [Hb), the transport 
spectrum now consists of the addition and removal levels 
U} = {4o,4h}, where ejg = E^ - = 6, and ejg = 
Eq — Eq = e + U. The expressions for the steady state 
currents [l5|,[20,[2^ through the hot contact (say) a = H, 
based on Egs. p^ and pT)) . are given by: 

J" = -,(i?o^iP"-i?f_,oTi+i?f,,pi-i?f_iT^), 
= (6-MH)(Po^i^^"-^f-.o^') 

where P^'s are the occupation probabilities of the many 
body state with 0, 1 or 2 electrons. We have dropped the 
index i within each N electron subspace, because only 
ground states exist within the framework of our spin- 
degenerate single orbital system. The solution for the 
set of master equations for this system based on Eq. (|13p 
is straightforward, and yields the following expressions 
for the occupation probabilities: 

P" = ^Pi^oP2^i, 
^ _ 1 

P = ^Po^lP2^1, 
2 1 

P" = ^Po^iPi^2, (28) 

with n being the normalization factor that ensures the 
sum of probabilities to be equal to unity. Here the total 
rate Ri^j = J2a-P?-^j^ given by the sum of the rates 
due to each contact a ~ H,C. To be specific, the addi- 
tion rates due to contact a ~ H in Eq. ([27|l are given by 
Rq^i = iHfie^-y'H), and 1h fif^ + U - ^ih)-. and 

the removal rates are given by Pfl^g = jh (l^/(e~Mif))- 
and Pf_,i = 7ff (1 -f{e + U^ fin)). 

We now plot the power (Eq. ipOl) ). and efflciency 
fEa. (|19p ) in the operating region in Fig. [31 In compar- 
ison with the non-interacting case, the domain of the 
operating region in the finite U case (green circles) is 
slightly different. This is because a finite U introduces 
a transport channel at e -I- C/ in addition to the already 
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FIG. 3: (Color online) Power and efficiency at finite U for 
■qc ~ 0.23. a) Power extracted in the operating region. The 
span of the operating region in the case of U — ksTc (green 
circles) can be different from that of the non interacting case 
(black solid). In general, the quantity Jq (see inset) is not 
identically zero when the electric current vanishes, b) Vari- 
ation of the efficiency in the operating region for different 
values of U: (i) U = (black solid), (ii) U — ksTc (green 
solid) and (iii) U = 2.5kBTc (gray dashed). Note that, with 
finite U such that U > hT, the efficiency is identically zero 
when the electric current vanishes under open circuit condi- 
tions {Vapp = Vs). The efficiency also reaches a maximum 
Vmax at finite power operation. 



existing one at e as shown in Fig. [TJc). For very small 
values of the interaction parameter U, specifically when 
U hT, higher order tunneling processes may become 
relevant. Such processes may only be captured by a per- 
turbative expansion beyond the second order in the tun- 
neling Hamiltonian [26| . 

We now plot the variation of the efficiency along the 
operating region Vapp in Fig. ^h) for different values of 
U. The trend of the variation of the efficiency with finite 
U is similar to what was noted in the non-interacting 
case. The efficiency reaches a maximum rjmax before be- 
coming zero. However, we also note from Fig. |3ljb) that 
the abruptness of this variation is less stark in compar- 
ison with the non-interacting case. In other words the 
maximum efficiency r]max occurs well within the domain 
of finite power. The introduction of interactions therefore 
also results in maximum efficiency within a finite power 
operation. This observed trend of the efficiency with ap- 
plied voltage as noted in Fig. El^b) may be qualitatively 



understood by analyzing the variation of currents with 
the appfied voltage Vapp. Based on Ea. (P7|) . one may 
recast an expression for the currents as: 

J = -«7(,/i(e) + J2(e + [/)), 
Jq = [e- fiH)Jiie) + ie + U - ^LH)J2{e + U), 

(29) 

where Ji(e) and J-zie + U) denote the contribution to the 
electric currents due to the transport channels at e and 
e + U and are given by: 

Jii^) = ^"^''^'"\ fH{e)-fc{e)), 
Me + U) ^ lin£^(^f^(^, + u)-fc{e + U)). 



(30) 



Likewise, 



(e - ^lc)Jlie) + {e + U- Mc) ^2(6 + U). (31) 



When Vapp is large enough to allow double occupancy in 
the quantum dot, the second transport channel e + U be- 
gins to conduct. The electrical current then redistributes 
between the two transport channels. From Eg. ([2^ . we 
note that the magnitude of Jq becomes more prominent 
as the contribution J2(e + U) increases. This causes the 
Jq to approach the zero value less rapidly with increasing 
bias than how the electric current would, thereby result- 
ing in an overall decrease in the ratio -ju- between them. 

-'q 

This causes the efficiency r] = ~ "^J^pp = JL to decrease 

9 " 

with increasing Vapp once the maximum r/max is reached. 
The applied bias at which this happens depends on U, 
and the above effect of the second transport channel will 
vary as U is increased. 

An important consequence of the introduction of this 
extra transport channel at e + U, is that both the cur- 
rents defined in Eq. ipT)) do not vanish at the same oper- 
ating point. Here, as shown in the inset of Fig. [3l[a) , Jq 
is finite even when the electric current J vanishes when 
{Vapp = Vs). The open circuit condition from Eg. ipn]) . 
can be deduced as Ji(e) ~ —Jzi^ + U). It then follows 
from Eqs.dSOl) and that ^ T^J^ = UJzie + U), 

and hence 



7^ 



J. 



H 



u 

7b 



U 



J2{e + U). 



(32) 



It is thus noted that similar to the non-interacting case, 
under open circuit conditions, a net entropy generation 
occurs, thus making spontaneous electron transfer pro- 
cesses irreversible. 

In order to further probe as to how the interaction 
U influences the achievable maximum efficiency, we plot 
the variation of rjmax (shown bold) with U in Fig. |4l 
We notice that with increasing U, the maximum effi- 
ciency reaches its global minimum around U ~ 2.7kBT, 
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FIG. 4: (Color online) Variation of the maximum efficiency 
with Coulomb interaction U. The maximum efficiency is equal 
to the Carnot efficiency for U — and asymptotically ap- 
proaches it when U ^ ksT. It reaches a minimum around 
U « 2.7fcsT. This variation is shown for a) rjc ~ 0.23 and b) 
r]c = 0.5. Also shown in each case is the comparison between 
the non-equilibrium calculation (bold) and that based on the 
figure of merit zT (brown dotted). Note that the difference 
between them becomes more prominent for larger values of 
rye or larger temperature gradients AT, thereby making the 
transport non-linear and hence the concept of zT less useful. 
The inset in (a) shows the variation of 1/zT with U for the 
chosen level configuration e — ~ 2kBTH at Vapp = 0. 



and asymptotically approaches Carnot efficiencies at very 
large values of U. As U is increased beyond U ~ 2.7kBT, 
the second transport channel e + U becomes less accessi- 
ble, and transport resembles the previous case with only 
one transport channel e. Thus, the important implica- 
tion here is that the variation of maximum efficiency 
with the introduction of interactions is non-trivial and 
non-monotonic. 

Our results are based on the evaluation of non- 
equilibrium currents and hence go beyond linear re- 
sponse. It is hence desirable to compare our results di- 
rectly with the conventional zT based evaluation which 
is valid only in the linear response limit. Using Eq.(I3|), 
the linear response maximum efficiencies calculated from 
zT are also plotted in Figs.jH^a) and (b), (brown circles). 
It must be noted from Fig. ^a) and (b) that the non- 
equilibrium calculation deviates from the zT based cal- 
culation [l^l (brown circles), and that this discrepancy 
is more pronounced for larger values of the Carnot effi- 



FIG. 5: (Color online) Variation of maximum efficiency with 
respect to 77c. a) The maximum efficiency (gray dotted) ap- 
proaches the Carnot efficiency and deviates more from the 
figure of merit zT based calculation (brown dotted) as the 
Carnot efficiency increases, b) Plot of the percentage de- 
viation of maximum efficiency between the non-equilibrium 
evaluation and the zT based evaluation. The maximum effi- 
ciencies at each value of ric here are taken from the respective 
global minimum {U ^ 2.7fcsT) in their variation with respect 
to U in Fig. H 



ciency. Also, comparing Fig. jH^a) and (b), it is seen that 
the deviation of the non-equilibrium efficiency from the 
Carnot efficiency with increasing U is less pronounced for 
larger values of the Carnot efficiency. 

To elucidate better, the discrepancy between the non- 
equilibrium evaluation and a zT based evaluation, we 
plot in Fig. [5j the variation of the non-equilibrium evalu- 
ation (gray squares) and the zT based evaluation (brown 
circles) of the maximum efficiency as a function of 77c, 
the Carnot efficiency. We note from Fig. ^a.) that the 
non-equilibrium calculation of rjmax deviates less from 
the Carnot value for both small and large values of rjc, 
with the maximum deviation in the intermediate re- 
gion. On the contrary, the zT based calculation devi- 
ates from both rjc and the non-equilibrium evaluation 
with increasing Carnot efficiency. From Eq.(I3]), in the 
zT based evaluation of the maximum efficiency, rjc is 
modulated by an increasing function of zT and is not 
strongly dependent on the operating conditions. Thus as 
r]c is increased, thereby increasing the applied tempera- 
ture gradient, non-equilibrium effects become prominent 
and transport cannot be adequately captured by the zT 
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FIG. 6: (Color online) Comparison between the efficiency at 
maximum power and various other limits for the interacting 
quantum dot set up. The efficiency at maximum power is eval- 
uated at f/ ~ 2.7kBT, which corresponds to the maximum de- 
viation from 77c in Fig. U The non-equilibrium evaluation for 
our set up using both symmetric electrostatic coupling (black 
diamond), and fully asymmetric electrostatic coupling (brown 
circles) is shown. The non-equilibrium evaluation assuming 
no electrostatic cou plin g to the cold contact, resembles the 
curves discussed in [IQI. The non-equilibrium evaluation of 
the efficiency at maximum power, in general, is dependent on 
the details of the set-up and need not be strictly bound by 
limits discussed in [131. 



based calculation. The inset in Fig. Hlja) shows the vari- 
ation of 1/zT with U, illustrating that zT — > 00 in the 
two opposite limits U ^ fc^T and U ^ fcsT. The per- 
centage deviation between the non-equilibrium calcula- 
tion and the zT based calculation of rjmax as a function 
of rjc is plotted in Fig.jSjb). 



C. Maximum power operation of an interacting 
quantum dot thermoelectric set up 

While Carnot efficiency poses the ultimate limit for 
any heat engine, there may or may not be other fun- 
damental limits involved under finite power operation. 
It has been shown that the maximum power operation 
of any Carnot engine is limited by the Curzon-Ahlborn 
efficiency ijca = 1 — VI ~ Vc [l^ ■ The study of the max- 
imum power operation of a non- interac ting quantum dot 
set up has been pursued previously [ll,0[l^. Here, we 
analyze the maximum power operation of the quantum 
dot system with the inclusion of Coulomb interaction. 

Analyzing maximum power operation implies maxi- 
mizing the electrical power output P = —JVapp. The effi- 
ciency at maximum power timp is calculated with respect 
to the operating point that maximizes the power output. 
We consider how this efficiency at maximum power t]]\jp 
compares with various limits discussed in earlier works 
(l6l lisl [l9| . Also, the efficiency at maximum power is 
evaluated for U w 2.7kBT, which corresponds to the 
maximum deviation of the efficiency from r]c in Fig. j?) 
We plot in Fig. [51 the calculation of the quantity r]M p 



under two set up conditions: (a) symmetric electrostatic 
coupling that is used throughout the paper (shown black 
diamond), and b) fully asymmetric set up in which the 
voltage applied across the cold contact is electrostatically 
decoupled to the quantum dot (shown brown circles) . In 
Fig. m we note that for smaller values of the tempera- 
ture difference, and hence smaller values of the Carnot 
efficiency ?7c, the efficiency at maximum power rjMP re- 
mains close to the Curzon Ahlborn limit and is approx- 
imately linear. In this limit, the curves follow a linear 
law. An important observation is that, similar to what 
was inferred in p^ . the efficiency at maximum power 
r]MP need not be bounded by the Curzon-Ahlborn effi- 
ciency for larger values of the Carnot efficiency rjc , and 
may indeed be larger. This questions the regime of ap- 
plicability of the Curzon-Ahlborn limit, which may only 
be valid for working conditions close to linear response. 
Consider an expansion for the Curzon-Ahlborn efficiency 
rjc A = 1 — 1/(1 — rjc) in powers of rjc = written as 



■nc_ 

2 



(33) 



from which it can be noted that for smaller values of 
AT and hence smaller values of the Carnot efficiency rye, 
the non-equilibrium efficiency follows the linear term af- 
ter which the quadratic term dominates. Notably, the 
deviation of the non-equilibrium efficiency at maximum 
power with respect to the Curzon-Ahlborn limit in Fig. [6] 
elucidates the fact that this limit need not be a funda- 
mental limit as the Carnot limit is. Physically, this im- 
plies that under non-equilibrium conditions, the leading 
term in the power expansion for "qcA deviates from a 
non-equilibrium evaluation, and importantly is specific 
to the set up. It has been pointed out in a recent work 
[l9| that in the limit of low dissipation, the efficiency at 
maximum power t^mp for a Carnot engine is bounded as 
< Vmp < V+, where -q^ = ^ with 77+ = (jz^- 
Note that our curve of the efficiency at maximum power 
is also not necessarily bound between the above two ex- 
trema. 

We thus note that the trend of the efficiency at maxi- 
mum power shown in Fig. |6]is similar to that of the non- 
interacting case analyzed in previous works [TgI . [TsI . [19| , 
when the quantum dot is electrostatically decoupled with 
the cold contact. The fact that the efficiency at max- 
imum power, under these conditions can approach the 
Carnot limit at certain larger values of AT (and hence 
77c), points out to the possibility of high power opera- 
tion at high efficiencies. While a large ratio of is not 
feasible at higher operating temperatures, it may be an 
interesting possibility in low temperature applications. 



IV. CONCLUSIONS 

In this paper, we analyzed the performance of an in- 
teracting quantum dot thermoelectric set up. This study 
was based on the evaluation of power and efficiency from 
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the non-equilibrium currents in the sequential tunneling 
limit. The operating region of the thermoelectric set up 
was identified and a general trend of the efiiciency in 
this operating region was identified. We showed that 
the much discussed aspect of reversible operation with 
Carnot efficiency, under open circuit conditions, in the 
case of non-interacting single orbital quantum dot sys- 
tem, only occurs in the limit of vanishing coupling to the 
contacts. In a general case, the efficiency reaches a max- 
imum in the operating region before dropping to zero at 
the open circuit operating point. In the non-interacting 
single orbital case, the efficiency can become very close 
to the Carnot value, if the coupling to the contacts is suf- 
ficiently weak. In the interacting case, we showed that 
this trend depends non-trivially on the interaction pa- 
rameter U. We also pointed out the clear discrepancy 
between our non-equilibrium evaluation of the maximum 
efficiency rjmax and the figure of merit zT based calcu- 
lation, which is only valid in the linear response limit. 
Comparison of the efficiency at maximum power with 
the Curzon-Ahlborn limit and other related bounds were 
also discussed. Here, it was shown that the inclusion 
of Coulomb interactions did not alter the already noted 
conclusions in the non-interacting case [l^, [l^ [l^ . How- 
ever, the trend of variation of the efficiency at maximum 
power is set up dependent. Our current theoretical treat- 
ment, however, is in the limit of weak coupling to the con- 
tacts, and symmetric contact coupling. In the regime of 
asymmetric and strong contact coupling, we expect novel 

tsics that may be introduced by asymmetric charging 
[35} to affect the thermoelectric transport processes. 
This will be an object of future research and possible ex- 
tension of the current work. 
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V. APPENDIX I: REVERSIBLE OPERATION 

The following important aspects of the so called re- 
versible operation in the case of a non-interacting quan- 
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